Reentrant Formation of Magnetic Polarons in Quantum Dots 
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We propose a model of magnetic polaron formation in semiconductor quantum dots doped with 
magnetic ions. A wetting layer serves as a reservoir of photo-generated holes, which can be trapped 
by the adjacent quantum dots. For certain hole densities, the temperature dependence of the 
magnetization induced by the trapped holes is reentrant: it disappears for some temperature range 
and reappears at higher temperatures. We demonstrate that this peculiar effect is not an artifact of 
the mean field approximation and persists after statistical spin fluctuations are accounted for. We 
predict fingerprints of reentrant magnetic polarons in photoluminescence spectra. 

PACS numbers: 75.50.Pp, 73. 21. La, 75.75.Lf, 85.75.-d 



Long spin memory timeji giant magnetoresistance^ 
robust magnetic ordering, and its versatile control in 



epitaxial and colloidal quantum dots (QDs) -• — are at- 
tributed to magnetic polarons (MPs), known for fifty 
years in bulk semiconductors The MP formation can 
be viewed as a "cloud" of localized spins, aligned through 
exchange interaction with a confined carrier spin. While 
the seminal studies of MPs in the bulkier— assumed that 
an impurity binds only one carrier, many experiments 
in QDs demonstrate multiple occupancies*^ - — We show 
that varying QD occupancy has important consequences 
for MP formation. It is conventionally understood that 
MPs form at low temperature (T) and vanish at high T, 
owing to thermal fluctuations of Mn spins. Here, we pro- 
pose an unexpected scenario where the temperature can 
enhance, rather than quench, MP formation and lead to 
reentrant magnetism i 20 ' 21 

We formulate a model of MPs formed in II- VI semicon- 
ductor QDs doped with Mn ions and apply it to simulate 
photoluminescence (PL) spectra^^^r— In a PL experi- 
ment^ the number of carriers captured by QDs depends 
on the total density of the photoexcited carriers, deter- 
mined by the laser intensity. Focusing on holes in type-II 
QDs^ we find nontrivial dependencies of MP binding en- 
ergies on both the laser intensity and T. 

Figure Q] shows our model. Epitaxial QDs typi- 
cally reside on a two-dimensional (2D) wetting layer 
(WL) in which electron-hole pairs are created by inter- 
band absorption of light. 26 Under illumination, a quasi- 
Fermi distribution is established in the WL 2D hole 
gas. Since the total area of the QDs is low relative 
to that of the WL, the quasi- Fermi level, fi(T,p) — 
+k B T\n (exp [pirh 2 / (m* h k B T)\ -l), is pinned by the 
non-equilibrium WL holes with 2D density p, where m* h 
is the heavy- hole effective mass. The population of the 
QDs by the captured holes is controlled by p. This as- 
sumption is justified by the fact that in many QD sys- 
tems, the times of capture and intra-dot relaxation are 
much shorter than carrier radiative recombination time. 
Therefore, quasi-equilibrium can be established in the 
valence band.— 

The exchange interaction of heavy holes with Mn spins 
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FIG. 1. (Color online) (a) A scheme of QDs grown on a wet- 
ting layer (WL), which is excited by light (cladding layer not 
shown). Type-II conduction/valence band (CB/VB) profile of 
a II- VI QD doped with Mn spins (green) . e and Ecv are the 
confinement and band gap energies. Hole quasi-Fermi level 
p.{T) lies in the continuum of WL states for T = 0. (c) QD 
states with corresponding energies, occupancies, and Mn-spin 
alignment: E\ state has zero energy and QD occupancy; £2,3 
state has one hole with spin +3/2 (—3/2), A is the exchange 
splitting; and E4 has two holes with a repulsive Coulomb en- 
ergy U. 



in a flat QD is highly anisotropic and described by an 
Ising Hamiltonian, 8 
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(1) 



where j3 is the exchange coupling constant, (Rj) and 
Szi (S Z j) are the position and the spin projection of the 
carrier (Mn)j2& The exchange interaction results in the 
Mn-magnetization, M z , and exchange splitting A of hole 
levels. We assume uniform hole wavefunction throughout 
the QD volume £1^2- thus we can relate M z and A = 
PM z /gpB- The maximum M z is M max = XMnNoflSgps, 
where XMn is the Mn fraction per cation, Nq is the density 
of cation sites, S = 5/2, g = 2.0 is the g factor, pb is the 
Bohr magneton, and Nq(3 ~ —1.0 eV^S 



The Gibbs free energy of the system is expressed as 



G sys (£) = G M „ (0 + ^(0, 



(2) 



in terms of the order parameter is £ = M z /M ma ^, where 
GW(£) is the Mn-spins contribution, and Fh(£) is the 
hole grand canonical (GC) free energy. Gmii(£) can be 
obtained by expressing the free energy of the Mn spins as 
a function of an external magnetic field and then applying 
a Legendre transformation ! 20 ! 30 



Gmu(0 = k B TN M -n 



In 



S mh[(l + l/2S)Bs\Q] 
sinh [Bg 1 (0/25] 
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where Bg 1 ^) is the inverse Brillouin function,— and 
A^Mn the number of Mn in the QD. Fh(£) is obtained 
from the QD states with 0, 1 or 2 holes [see Fig.pTc)]. For 
transparency, we consider only one nonmagnetic single- 
hole level and neglect the possibility of magnetization in 
the presence of two holes This yields 
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where e is the single-hole confining energy, U is the re- 
pulsive Coulomb (charging) energy, and the maximum 
splitting is A max = XMn\No/3\S. To elucidate some 
interesting phenomena, we first use a standard mean 
field (MF) approximationji^— commonly used by many 
authors By numerically minimizing Eq. @, we 
obtain the value £mf- The MP energy, Emp, is de- 
fined as the expectation value of Eq. ((T|). We calcu- 
late it from E M p = -fc s TA max c? (lnZ sys ) /dA maxJ 2£ The 
MF partition function Z sys (£mf) = Z Mn (£mf) Z h (£mf) 
is expressed in terms of the Mn and hole contri- 
butions: Z M „(£mf) = e -GM„(?MF)/fc B T Z h (^MF) = 



,-, A max ^MF —(F—lA/kaT ■ ^ 

Emp = — „ tt r e (e sinh 
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In Fig. [2] MF predicts multiple phase transitions and 
reentrant magnetism. For low- hole densities (solid line), 
the system exhibits a second-order phase transition at 
Tc3 — 27 K. In Fig. [U(b), the location of minimum 
Gsys(C) continuously goes to zero, a signature of sec- 
ond order transition. This is similar to the usual MP 
case (single particle), since for any T < 30 K, the prob- 
ability of finding a single hole in the QD is dominant. 
For high- hole densities (dotted line), MF predicts a first- 
order transition at Tci = 5 K, consistent with the dis- 
continuous shift of the G sys (£) minimum to £ = at 
Tci in Fig. [2£c). At T < Tci, magnetism is present 
since the QD is occupied by 1, rather than 2 holes 
(despite /i ~ 1 meV in the continuum), according to 
E-2 — (i < E4 — 2/i [Fig. [He)]. This inequality is satis- 
fied as long as the ordering of Mn-spins sufficiently lowers 




FIG. 2. (Color online) (a) Reentrant behavior in the MP 
energy, seeE-MP [Eq. ([5J], and various critical temperatures. 
Free energy evolution reveals the second- (b) and first-order 
(c) phase transitions for p — 4x 10 9 and p = 9 x 10 10 cm" 2 , re- 
spectively. Parameters, partially guided by (Zn,Mn)Te/ZnSe 
QDs4 e= -36 meV, U = 30 meVr^ A max = 69 meV (see 
Fig. HJl, x Mn =2.6%, No/3 = -1.05 eV, and m* h = 0.21. QD 
volume, fl, is cylindrical with a 50 A radius and a 25 A height. 



the energy of the single hole state E^. When T > Tci, 
E-2 — [a > E4 — 2/i and the QD becomes doubly occupied 
suppressing MP formation. Above Tci, magnetism does 
not reappear since thermal Mn excitations completely 
quench magnetic order, before /i approaches e + U to 
promote single occupancy. 

The QD exhibits reentrant magnetism at p — 5 x 
10 cm" - 2 (dashed line). For T < Tg 2 , the scenario is 
the same as the dotted line (X^ plays the role of Tci)- 
At T > Tg 2 , with E2- M < E 4 - 2/i, the QD is singly 
occupied as a result of fx moving quickly towards e, thus 
MP reappears by a first-order transition. The transition 



at T = Tp 3 has the same origin as the solid curve. For 
T 3> Tq 3 , the dot becomes emptied. 

However, MF theory neglects the possibility for the 
system to deviate from the equilibrium value, £mf- This 
leads to unphysical thermodynamic phase transitions in 
small systems. Could this also imply that the described 
reentrant magnetism is only an artifact of the MF the- 
ory? To address this question and better understand the 
validity of the behavior predicted at the MF level, we 
formulate a fluctuation approach (FA). Statistical fluc- 
tuations are included in the partition function by inte- 
grating over all possible values of the order parameter. 12 
Correspondingly, we employ Z sys = f_ t e~ G 'y^/ kBj d£ 
to implement the framework used for MF [recall Eq. (JSJ], 
and obtain the average exchange energy 

A m 
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This Emp for the GC ensemble is similar to that of the 
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canonical ensemble^ However, Z sys now contains multi- 
ple occupancies and the numerator is weighted by e >i ' kBT ', 
which decreases with increasing T [Fig. H£b)]. 

We are now able to directly compare MF and FA re- 
sults. The sharp MF phase transitions [-Emp hi Fig. [5], 
become smeared out, as seen in Fig. [31 FA yields finite 
Emp at any finite T. This is expected from averaging 
of G sys over £, implicit in Eq. ([B]), including strong con- 
tributions from the competing local minima [Fig. [2jc) 
at Tci]- For example, MF reentrant magnetism from 
Figs. [Ua) and[3jb) for p = 5 x 10 9 cm~ 2 is absent in FA 
(see Fig. |3j, since the local minima of G sys at |£| > 
contribute strongly in the temperature range of 9—13 
K. Surprisingly, the dotted and solid curves in Fig. [3] 
show FA reentrant Emp even at room temperature, while 
for the same p no MF reentrant behavior was seen. In 
FA, the increase in |-Emp| at higher T is due to the the 
QD occupancy decreasing from 2 to 1 holes. The maxi- 
mum reentrant |Emp(T 1 )| coincides with the average oc- 
cupancy of 1. Even though inclusion of statistical fluctu- 
ations removes reentrance for some hole densities, it also 
yields a smoothed version of the same effect for higher T, 
where MF predicts Emp = 0. A peculiar nonmonotonic 
Emp(T) is therefore not limited to the MF description. 

Furthermore, the reentrant MP is not restricted to the 
above parameters, but occurs for a range of iMn, £, U, 
and p. For example, we find reentrant Emp for 1.5 x 
10 10 cm -2 < p < 2 x 10 12 cm~ 2 , with other parameters 
fixed. Conversely, the reentrant MP is present for 1.5% < 
iMn < 2.6%, if the remaining parameters are fixed. 



energy ~ Eqv > C 36 ( see Fig.[T]) and with opposite spins. 
The total spectrum, J tot is the superposition of the lines 
generated by the 2->l, and 1 — » transitions, 



p=5x10 cm 
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FIG. 3. (Color online) (a) Reentrant magnetism in the 
fluctuation approach (FA) [see Eq. fl6}] forp = 5xl0 9 , 9xl0 10 , 
4 x 10 11 cm -2 , (b) Comparison of the FA and mean field 
(MF). FA (dot-dashed line) shows no reentrant magnetism 
predicted for MF and p = 5x 10 9 cm -2 (dashed line). Other 
QD parameters are from Fig. [21 



We next discuss how reentrant magnetism could be 
observed in PL experiments. To obtain the PL spectrum, 
we assume low QD density and sufficiently high intensity 
of exciting light for two CB electrons to always be in 
the vicinity of each QD, occupying the lowest possible 



-ftotM = /i-kjM +/2->-i(u>)- 



(7) 



We assume that the Mn-configuration does not change 
during a recombination event. The intensity of each PL 
line is / = J2i fPi w ifd ~ {E% — Ef)], where Wif is the 
transition rate, pi is the thermodynamic probability of 
the initial state, fku is the energy of the emitted photon, 
and Ef (Ei) is the energy of the final (initial) state of 
the system. We replace the above with J d£. For 1 
— > transitions, the system is in an initial state with a 
hole of spin up (down), which later recombines with a 
spin-down (up) electron^ 2 . The intensity of this line is 



h^o(X) = c (T)e-^-^lkBT e -x/k B T 

x e -G Mn (2X/A nlax) /k B T e (Amax _ | 2X | 



(8) 



Here, X(u>) = fvjj — Eqv — s is the shifted frequency, 
c(T) oc (f2//3)|<i ct ,| 2 /iJ S y S , d cv the dipole matrix element, 
and 8 is the step function. For the 2 — > 1 transition, 
there are two holes and two electrons of opposite spin in 
the initial state. The intensity is 

l2-n(X) = c(T)e- {2e+u - 2 ^ /kBT e~ GM » {2{x - u)/A ™ x)/kBT 
x9(A max -2\X-U\). (9) 

In Fig. [4l the PL spectrum shows the evolution of the 
peaks for transitions 1 — > 0, centered at A < 0, and 
2—^1, centered at the charging energy, X = U, since 
hw = Ei - Ef = (2e + U + E cv ) - e , 37 From Eqs. ©, 
©, and ([9]) it follows that the Mn-contribution to the PL 
is T independent: the amplitude of the 1 — > (2 — > 1) 
peak at different T is proportional to the probability of 
finding a single (double) occupied QD. For T < 5 K, 
the 1 — > line dominates, while it becomes negligible at 
~ 10 K where the system is virtually nonmagnetic and 
the double occupied state (2 — > 1 line) is dominant. At 
higher T, due to the shift of ft toward CB, the probability 
of single occupancy increases, and for T > 100 K, the 
probability of zero occupancy increases. The resulting 
T-dependencies are remarkably nonmonotonic for both 
/tot peak position (red and blue shifts), and I to t peak 
intensity; a signature of reentrant MPs, consistent with 
£mp(T) in Fig. 01 

What are the semiconductor systems where the reen- 
trant MP could be found? Recently, a nonmonotonic PL 
red shift was observed in type-II (Zn,Mn)Te/ZnSe QDs, 4 
which have partially guided our choice of parameters. 
However, the reentrant magnetism should not be lim- 
ited to type-II systems. The necessary condition is the 
T-dependent multiple occupancy, readily seen in both 
type-I I 16 i 17 and type-I QDsi 15 ' 38 Multiple occupancy can 
be reached by raising photo-excitation intensity, which 
may first lead to weakening of MPs (blue shift) through 
Mn-spin heating] 22 ! 39 " — Nevertheless, an increased blue 
shift was attributed to double occupancy in type-I mag- 
netic QDs^ 
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FIG. 4. (Color online) (a) A PL spectrum for reentrant mag- 
netism showing 2 — > 1 and 1 — > hole occupancy transitions 
(p — 9 x 10 10 cm -2 ). The thick lines emphasize important 
features in the PL spectrum, while the thin lines between 
them show their evolution at intermediate equidistant T. (b) 
Overall PL peak position marked with filled circles. As T in- 
creases, 1 — > peak shifts to zero energy. QD parameters are 
from Fig. [2] 



Considering only a reduced T-range could conceal the 
presence of reentrant MP. An initial steep decline in 
£mp(T) (see Fig. [3]) is similar to conventional MPs, while 
a slightly higher T region (e.g., 5 to 30 K for solid line in 
Fig. [3]) could be misinterpreted as a final thermal breakup 
of MPs. Thus, further experimental studies of power and 
T dependence are important. Specifically, it would be de- 
sirable to consider single self- assembled QDs&^2 to reduce 
uncertainties due to inhomogeneous averaging, and focus 
on moderate xmh to suppress the Mn-Mn antiferromag- 
netic interactions.- Colloidal (II,Mn)VI QDs&i£ showing 
a robust MP formation are also promising candidates to 
test some of our predictions. 

Even though in this work we have only focused on 
reentrant magnetism, we expect that further studies of 
the implications of multiple occupancies will lead to ad- 
ditional surprises in both epitaxial and colloidal QDs. 
Since prior findings in QDs were successfully applied to 
different finite fermion systems^ it may be possible to 
seek other promising paths for observing reentrant mag- 
netism. 
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